Mixed state Pauli channel parameter estimation 
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The accuracy of any physical scheme used to estimate the parameter describing the strength of 
a single qubit Pauli channel can be quantified using standard techniques from quantum estimation 
theory. It is known that the optimal estimation scheme, with m channel invocations, uses initial 
states for the systems which are pure and unentangled and provides an uncertainty of 0(l/s/m). 
This protocol is analogous to a classical repetition and averaging scheme. We consider estimation 
schemes where the initial states available are not pure and compare a protocol involving quantum 
correlated states to independent state protocols analogous to classical repetition schemes. We show, 
that unlike the pure state case, the quantum correlated state protocol can yield greater estimation 
accuracy than any independent state protocol. We show that these gains persist even when the 
system states are separable and, in some cases, when quantum discord is absent after channel invo- 
cation. We describe the relevance of these protocols to nuclear magnetic resonance measurements. 

PACS numbers: 03.65.Ta, 03.67.-a,03.65.Ud 



I. INTRODUCTION 

Quantum estimation considers situations in which a 
quantum system or evolution operation depends on one 
or more parameters, the values of which are to be es- 
timated as accurately as possible. Such estimation will 
have to be accomplished by subjecting physical systems 
to physical procedures and consequently the accuracy of 
estimates will be governed by the relevant laws of physics. 
This realization has resulted in substantial study of con- 
straints and quantum advantages available for estimation 
in the context of quantum theory and has yielded results 
with important implications for quantum metrology [1- 
13]. 

In general, any quantum estimation procedure consid- 
ers a quantum operation of known form but which de- 
pends on one or more parameters. The focus of this arti- 
cle is the Pauli operation (or channel), T(A), which acts 
on a single qubit. This transforms a system, initially in 
the state p, via 



f (A) 

p H> pi (A) := (1 - A) p + Xa n pa n 



(1) 



where ^ A ^ 1 and a n is a Pauli operator with n repre- 
sent a direction in space. For example, n = x represents 
the bit-flip operation and n = y is the phase-flip oper- 
ation. The assumption is that the direction n is known 
but the parameter A is unknown. The task is to esti- 
mate A with minimal uses of physical resources, such as 
qubits, number of applications of the channel or ancillary 
quantum operations. 

Channels of this type are important for the following 
reasons. First, the three channels for which n = x,y and 
z are sufficient to generate the types of errors that are 



considered in quantum information processing [14]. Esti- 
mating their strengths via the parameters will be relevant 
for practical quantum information processing protocols. 

Second, specific examples of such channels appear in 
contexts such as nuclear magnetic resonance (NMR), in 
which the phase-flip operation generates the dephasing 
(or T 2 ) processes [15, 16]. Here, the density operator of 
a single spin- 1/2 system evolves with time t as 
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(2) 



where T 2 is a constant with units of time. This is related 
to the phase-flip parameter by 



A=(l-e 



-t/T 2 



)/2. 



(3) 
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Measuring or estimating the phase-flip parameter is thus 
equivalent to measuring T 2 . 

Third, quantum parameter estimation can illustrate 
fundamental differences between classical and quantum 
worlds. For example, consider two methods of estimating 
the parameter governing unitary phase shifts [8] . Apply- 
ing the unitary phase operation once to each of n inde- 
pendent quantum systems, yields an optimal uncertainty 
of 0(l/y/n), also called the standard quantum limit. 
This is analogous to a classical repetition and averaging 
procedure. A method which applies the unitary phase 
operation once to each of n quantum systems which are 
initially in a particular entangled state yields an uncer- 
tainty of 0(l/n). This is the Hciscnberg quantum limit 
and is optimal for this type of operation [8, 17, 18]. Ev- 
idently the resource of entanglement, occurring in quan- 
tum systems, has yielded an advantage here. 

Parameter estimation for Pauli channels has been con- 
sidered in the context of finding an optimal scheme us- 
ing all possible input states of the quantum systems in- 
volved [7, 19-23]. The results, reviewed in section III, 
are that the optimal accuracy can be attained when the 
channel is invoked multiple times on independent and 
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uncntanglcd systems, each of which is initially in a pure 
state [19]. The manifestly quantum resource of entangle- 
ment has no role to play in this circumstance. However, 
in situations such as solution state NMR, initial states 
are highly mixed. The question remains of whether an 
independent channel invocation approach is still optimal 
when the only initial states available are mixed. This 
variant of quantum parameter estimation has been in- 
vestigated for unitary operations [24-27] where advan- 
tages of the type using entangled pure states survived 
with mixed, correlated initial states. 

In this article we consider the question of Pauli channel 
parameter estimation when the initial states of the avail- 
able quantum systems are mixed. We aim to compare the 
accuracy of an independent channel use approach to one 
where correlated or entangled states are used. This will 
be important for Pauli channels in solution state NMR 
and other situations where the initial states are inher- 
ently mixed. Particularly, for NMR it is not practical 
to purify the state by a process of measurement followed 
by selection based on measurement outcomes; the full 
range of available initial does not include pure states. 
Our approach is similar to that of Modi, et. al. [26] but 
the results arc subtly different, particularly in terms of 
the relationship between the presence of entanglement or 
discord when improved accuracy is attained. 

We shall focus on the phase- flip channel but the results 
hold for any Pauli channel. Section II reviews the gen- 
eral framework for quantum parameter estimation and 
general bounds for phase-flip parameter estimation. Sec- 
tion III describes initial states available and gives the 
optimal independent channel use estimation accuracy; 
this provides a "classical" benchmark against which other 
schemes can be compared. Section IV outlines a partic- 
ular approach that uses correlated quantum states and 
gives central new results which form the foundation our 
the rest of the article. Sections V and VI describe two 
cases of the main result with notably different outcomes. 
Finally section VII discusses the presence of correlations 
such as entanglement or quantum discord in our estima- 
tion schemes. 



II. QUANTUM PARAMETER ESTIMATION 

Typically a physical quantum estimation procedure 
commences with various quantum systems prepared in 
known initial states. The parameter dependent operation 
is applied multiple times to some of these systems, pos- 
sibly interspersed with ancillary parameter-independent 
quantum operations. The physical procedure terminates 
with measurements on the individual quantum systems. 
The parameter is inferred by data processing from the 
measurement outcomes. Fig. 1 illustrates a generic quan- 
tum estimation scheme which uses m invocations of the 
operation and n quantum systems. The general proba- 
bilistic nature of the outcomes of measurements on quan- 
tum systems means that the parameter inferences will 




FIG. 1. Schematic of a quantum estimation procedure using n 
quantum systems, each represented by a horizontal line. This 
procedure uses m invocations of the parameter-dependent op- 
eration F(A), interspersed with parameter-independent uni- 
tary transformations Ti, . . . , T m . The initial state of all quan- 
tum systems is pi and the pre-measurement state /5f(A). Mea- 
surements on individual systems follow the final unitary T m . 



fluctuate about the true parameter value. A key goal of 
any parameter estimation scheme is to minimize these 
fluctuations. 

Various estimation schemes use differing resources (e.g. 
number of qubits, number of channel invocations, num- 
ber of two qubit entangling unitaries) and a fair compar- 
ison requires precise specification of the resource costs. 
A common assessment of the resource cost uses a black- 
box query count within a quantum circuit [8, 17, 18, 28] 
such as illustrated in Fig. 1. In this approach, used in 
this article, the only relevant resource cost is the num- 
ber of times, to, that the phase-flip operation appears 
in the quantum circuit for the protocol. We will only 
compare schemes that use the same number of channel 
invocations. 

In the scheme of Fig. 1, p\ can be assumed to be a 
product of states for the individual qubits since an po- 
tential entangled state can be generated from this via 
an appropriate unitary, which can be subsumed into To. 
Given this, estimation schemes can broadly be grouped 
into two classes. 

In the first, called an independent channel use protocol, 
the channel is applied once to each of to qubits, which 
arc initially in a product state. This is followed by mea- 
surements. In the context of Fig. 1 this means that all of 
{Ti, . . . , T m } are non-entangling amongst the qubits (a 
non-entangling unitary transforms any state which is a 
product of states for individual quantum systems to an- 
other such product state). In this protocol To is a prod- 
uct of single qubit unitaries while T± , . . . , T m are products 
of single qubit unitaries sandwiched between two qubit 
swap gates. Throughout the procedure of this proto- 
col the state of the entire system is always a product of 
states for the individual qubits and there are no quan- 
tum correlations present. The independent channel use 
protocol is analogous to the classical strategy of invoking 
the physical operation independently on multiple physi- 
cal systems; this generally enhances estimation accuracy. 

In the second approach, called a correlated state proto- 
col and only available to quantum systems, at least one 
of {Ti,...,T m } is an entangling unitary and, at some 
point during the procedure, the system will not be in a 
product state. 
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The task of schemes falling into either class is to find 
admissible unitaries {Ti, . . . , T m } that maximize the esti- 
mation for fixed resource costs. The question is whether 
there is some correlated state protocol that provides bet- 
ter estimation accuracy that any independent state pro- 
tocol. If so, it would indicate that some quantum re- 
sources, possibly generated by an entangling unitary, 
have been harnessed so as to outperform any estima- 
tion protocol merely analogous to the classical strategy 
of repetition, as employed by an independent channel use 
protocol. 



A. Quantum Fisher information 

Quantifying the accuracy of quantum estimation pro- 
cedures is done using the following general scheme. De- 
noting the outcome of the measurement on the j th system 
by Xj, the collective outcomes of the physical procedure 
of Fig. 1 can be represented by {x\, . . . , x n }. The param- 
eter is estimated by applying these an estimator function, 
A C st, resulting in the estimate, A = A cs t(xi, . . . , x n ). As 
a result of statistical fluctuations in measurement out- 
comes, inherent in quantum systems, the same estimator 
will yield possibly different results from one run to an- 
other which uses an identical sequence of quantum op- 
erations. The statistical accuracy of the estimate will 
therefore rely on the choice of estimator together with the 
probability distribution for the measurement outcomes. 
One means of quantifying the accuracy of this estimate 
uses the mean square error, 



mse (A) := ( I A — A 



(4) 



where the angle brackets indicate the mean over all pos- 
sible measurement outcomes. The possible estimators 
consider in this article will be unbiased estimators, i.e. 
(A) = A, and here mse (A) = var(A), where the variance 
in the estimate is 



var (A) 



A -(A) 



(5) 



The variance is the fundamental entity which quanti- 
fies the accuracy of the estimate. A fundamental re- 
sult [9, 11, 29] is that, regardless of the estimator func- 
tion, the variance is bounded from below via the Cramer- 
Rao bound 



var (A) > 



1 



F(X) 



where F(X) 
F(X) := 



d]np(xi, . . .x n \X) 



n 2 



dX 



(6) 



dxi...dx n . (7) 



is the Fisher information associated with the probability 
distribution (conditional upon A) for the measurement 



outcomes, p(x\, . . . x n \X). If each measurement outcome 
is one of a set of discrete possibilities, the integral is re- 
placed by a summation. 

It can be shown that under general conditions, there is 
always an estimator [29] , which asymptotically saturates 
the bound of Eq. (6). Thus in classical parameter esti- 
mation the Fisher information quantifies accuracy of the 
estimate and when comparing estimation schemes, that 
with the greater Fisher information is more successful 
than any competitor with a lower Fisher information. 

Given a particular estimation procedure with Fisher 
information, F(X), an obvious and conventional classical 
strategy for improvement is to repeat the identical proce- 
dure multiple times independently. Mathematically this 
is equivalent to sampling the same probability distribu- 
tion in an independent and identical fashion. If this is 
done to times then the Fisher information is mF(X); this 
governs the typical behavior as a function of m in classical 
parameter estimation [11]. For example, in the context 
of parameter estimation in classical physical systems, the 
same physical procedure can be applied to to identical 
systems. The resulting mean square error is reduced by 
a factor of to. 

For parameters pertaining to the evolution of quantum 
systems, the probability of measurement outcomes is de- 
termined by: (i) the choice of measurement and (ii) the 
choice of initial state pi, which determines the final state 
Pf(A). Therefore, the classical Fisher information for any 
given final state depends on the choice of the measure- 
ment. However, considering all possible quantum mea- 
surement procedures on the system and any possible an- 
cillary quantum systems shows that the classical Fisher 
information satisfies the quantum Cramer-Rao bound, 



F(X) ^ H(X) 



(8) 



where the quantum Fisher information H(X) is indepen- 
dent of the choice of measurement [5, 9, 11]. The quan- 
tum Fisher information depends only on /5f (A) and is cal- 
culated via [9] 



H(X)=Tr Pi (X)L 2 (X) 



(9) 



where the symmetric logarithmic derivative (SLD) or 
score operator, £(A) is defined implicitly via 



dX 



L(A)p f (A)+pf(A)L(A) 



(10) 



A useful alternative to Eq. (9) is 
H(X)=Tr\ Pi (X)L 2 (X) 



2 



= Tr 



L(X)p i (X)L(X)+p i (X)L 2 (X) 
dp { (X) 



OX 



-L(X) 



(11) 



There exist methods [9] for computing the SLD based on 
eigenvalue and eigenstate decompositions of pt(X) and 
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conditions for the existence of measurements which sat- 
urate the bound of Eq. (8). In general there does ex- 
ist a measurement procedure which asymptotically at- 
tains this bound [30] . Quantum Fisher information then 
serves as a measure of the possible accuracy of the es- 
timation procedure. The actual measurement and data 
processing of the outcomes are irrelevant here and all that 
matters is the final state of the system, /5f(A). The cen- 
tral task in any quantum parameter estimation scheme 
is to determine the initial state for the quantum system 
which maximizes the resulting quantum Fisher informa- 
tion and, if possible, the associated measurement which 
yields a classical Fisher information equal to the maxi- 
mum quantum Fisher information. This reduces the task 
of designing quantum parameter estimation schemes to 
that of tailoring pt(X) so as to maximize the quantum 
Fisher information. 



B. Determining the SLD 

Computing the SLD can be done by computing eigen- 
states and eigenvalues of pt(X) but this generally diffi- 
cult. However, for certain operators, there are methods 
for computing the SLD which do not require the eigen- 
values and eigenstates. 

Proposition 1: Let A be a diagonalizable operator on 
a two dimensional vector space such that Tr A ^ 0. Let 
a := Tr (i 2 ) - (Tri) 2 . Then §| = (LA + AL)/2 where 



L = 



TrA 



dA <91n(Trvl) ~ 
2 dX d~X A 



(12) 



if a = and 



L =- 



1 



TrA 
d_ 

dX 



dA dlna ~ 
~dX dX~ 



In a - In (TrA) 



(13) 



if a 0. 

Proof: Denote the eigenvalues of A by a\ and a^- The 
characteristic equation for A is (A — a\I)(A — a^l) = 0. 
Expanding this gives 

A 2 - Tr [i]i + axa 2 I = 0. 



Further algebra yields 



.4 
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A* 



1 



al 



TrA V 2 
where we have used the fact that Tr A ^ 0. Thus 



(14) 



dA 
~dX 



1 9 In (Tri) ln 



Tri dX 
1 / dA 



(*-?') 



Tr A 



dx A + A dx ~ 2 dx , 



(15) 



If a = 0, Eq. (15) reduces to 
dA 1 din (Tri) 



dX 



TrA 



dX 



,4 J 



TrA 



dA j -dA 

dx A + A dx 



and by inspection, Eq. (12) holds. On the other hand, if 
a^0 then Eq. (14) can be inverted to give / = 2[i 2 — 
Tr (A) A] /a and Eq. (15) becomes 



d\n(TrA) 




TrA 



A— -lOHf 
~dX~2~dX 



A z - Tr (A) A 



(16) 



Inspection of Eq (16) shows that the SLD is given by 
Eq. (13). □ 

Some situations require SLDs for operators that act 
on vector spaces of dimension greater than 2. If such an 
operator can be decomposed into the sum of operators 
acting on orthogonal spaces then the quantum Fisher in- 
formation can be determined in a piecewise fashion. 

Proposition 2: Suppose that p — pk where the 
supports of fjj and pk are orthogonal whenever j ^ k. Let 
Ik be the SLD for p k , i.e. ^ = (L k p k + p k L k )/2. Then 
the SLD for p is L = ^ fc L k and the quantum Fisher 
information is H — J2 k Hk where H k = Tr (^-L k ). 

Proof Suppose that ^ = (L k p k + p k L k )/2. Then 



dp_ _ dp k 
dX~ ^ dX 

k 



L k p k + PkLk 



(17) 



The support of L k is either the same as or a subspace of 
the support of p k . Thus if j ^ k then Lj p k = and 



dp 
dX 



Lp + pL) 



(18) 



The result regarding the quantum Fisher information fol- 
lows by applying similar reasoning to Eq. (13). □ 



C. Bounds on the quantum Fisher information for 
phase-flip parameter estimation 

A useful tool for addressing bounds on the quantum 
Fisher information for the phase-flip operation, already 
considered elsewhere [7, 19, 21], involves a representa- 
tion of the channel using an additional ancilla qubit as 
illustrated in Fig. 2. 
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Po 



V(X) 



discard 



Pfs(A) 



FIG. 2. Unitary representation of a phase-flip channel. The 
lower qubit is the channel qubit and the upper an ancilla 
qubit, which is required to be in the initial state, |0). Af- 
ter evolution under the two indicated unitary operations, the 
ancilla qubit is ignored or discarded. 



This pair of qubits is subjected to the indicated uni- 
tary operations where the additional 11 coin-toss" unitary 
operation is 



V(X) := 



a/A 
a/1 - A 



-y/T= A 

Vx 



(19) 



Following the controlled- Z unitary, the ancilla qubit is 
discarded. Mathematically, this amounts to taking the 
partial trace over the state space for the ancilla qubit. 
It is straightforward to show that the resulting transfor- 
mation of the channel qubit is identical to the phase-flip 
operation. 

The quantum Fisher information of the phase-flip 
channel alone is bounded from above by the quantum 
Fisher information of the entire extended channel of 
Fig. 2 since the set of all possible measurements on both 
qubits includes the set of all measurements made on 
the channel qubit alone. When evaluating the quan- 
tum Fisher information of the extended channel the 
controlled- Z unitary and the initial state of the channel 
qubit can be collectively absorbed into a general measure- 
ment procedure since they are parameter-independent. 
The only relevant entities are the coin-toss unitary and 
the initial state of the ancilla qubit. Thus the quantum 
Fisher information for the extended channel is exactly 
that for the coin-toss unitary applied to the ancilla in the 
state |0). Proposition 1 and Eq. (11) apply here and give 
that the quantum Fisher information for the extended 
channel is 1/A(1 — A). 

Multiple uses of the phase-flip operation can be repre- 
sented by a similar extension involving multiple coin-toss 
unitaries and ancilla qubits, each in the initial state |0). 
According to the scheme above, these parts of the proce- 
dure occur independently of each other and it is straight- 
forward to show that the quantum Fisher information of 
the extended channel (and thus the upper bound on any 
estimation procedure using any channel input states) us- 
ing m phase-flip operations is m/A(l — A). This result is 
independent of the state of the channel qubits and also of 
any unitaries T\ , . . . , T m that are applied to the channel 
qubits between channel invocations. Thus for any param- 
eter estimation scheme involving m channel invocations 
the quantum Fisher information satisfies 



A(l-A)' 



(20) 



This absolute upper bound has appeared elsewhere in 
more general context [19]. 



III. INPUT STATE CHOICES AND 
CONSTRAINTS 

The quantum Fisher information partly depends on 
the choice of the initial state for the system. It will be 
essential to provide a general mathematical description 
of the various initial states; the quantum Fisher infor- 
mation and possible constraints on initial states will be 
expressed partly in terms of parameters appearing in this 
description. 



A. General product input states 

In general the initial state of the system can be as- 
sumed to be a product of states for individual qubits; 
initial entangled states can then be accommodated via 
appropriate choice of T . The general density operator 
for an individual qubit has form 



P := 



I 



(21) 
,) is a 



where r • a = r x a x + r y & y + r z a z and r = (r x 
conventional three dimensional real vector whose norm, 

or polarization, satisfies r = |r| := ^r x + r y + r' 2 z ^ 1. 

We assume that all qubits have the same initial polar- 
ization. The only remaining possible differences between 
initial states arises from the orientations of r for each 
qubit. By including appropriate single qubit rotations in 
To, any differences in these orientations can be removed. 
Thus we assume that r is identical for each qubit and 
therefore p- x = po <g> po ® ■ ■ ■ ® po {n factors) where 



Po 



I (;+,,) 



(22) 



for a particular given r. The freedom to include single 
qubit rotations in T thus implies that the only relevant 
input state parameter for the quantum Fisher informa- 
tion is the polarization, r. 



B. Pure input states 

A subset of possible input states are those which are 
pure, i.e. r = 1. Consider a single application of the chan- 
nel to a qubit in the initial state \+y) — (|0) + i\l)) 
with all Tj = I. Straightforward application of proposi- 
tion 1 and Eq. (11) yields a single use pure state quantum 
Fisher information 



H. 



s pure opt 



(A) 



1 



A(l-A)' 



(23) 
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identical to that of the bound of Eq (20) with m — 1 and 
therefore optimal. In fact, the projective measurement in 
the basis {\+y), \-y)} , where \±y) = (|0) ± i|l>) fy/2, is 
readily shown to yield a classical Fisher information equal 
to H s p Ure opt (A). 

It follows that if this is repeated independently with 
m channel invocations, each applied to a single qubit as 
described above, the resulting classical Fisher informa- 
tion is m times that for a single channel use. But this 
saturates the bound of Eq. (20) and thus the optimal 
pure state quantum Fisher information with m channel 
invocations is 



pure 



opt (A) — 



A(l-A) 



(24) 



Clearly for pure input states the optimal quantum Fisher 
information can be attained using product states (this is 
an example of the independent channel use protocol) . 
For pure input states no correlated state protocol can yield 
a better estimate than the particular independent chan- 
nel use protocol describe in this section. Evidently entan- 
glement or quantum correlations between qubits cannot 
yield any advantages for phase-flip parameter estimation 
with initial pure states. 



C. Independent Channel Use Protocol 

In the more general independent channel use protocol 
the initial states are not pure. Suppose that the operation 
is applied once to each of m qubits, each of which is in the 
identical initial state, given by Eq. (22) and which appear 
leftmost in the mathematical representation of the entire 
system state. The final state of the collection of qubits 
is pf = pi s ® • ■ ■ <8 Pf 8 ® p ® • ■ ■ ® p (to factors of pf s ) 
where the final state for the individual qubits to which 
the operation is applied is pf s = (1 — A) po + \o z p$<j z . 
Let L s be the SLD corresponding to pf s and let H S (X) be 
the resulting quantum Fisher information for this single 
qubit. It is straightforward to show that the SLD for all 

n qubits is L s ®i®---®i + i®L s ®i®---®I^ 

with one term for each qubit to which the operation is 
applied. Then Eq. (11) implies that the quantum Fisher 
information for the entire system is 



H{\) = m H S (X). 



(25) 



The quantum Fisher information for a single channel 
use on one qubit can be computed by first determining 
the SLD. The conditions of proposition 1 apply and, since 
Pf s is not a pure state, a = Tr (p 2 s ) — (Tr pf s ) ^ 0. Thus 



is 



Pi: 



Tr pf s pf. 



l-Tr[p f 2 J 



Pit 



(26) 



Substituting into Eq. (11) gives 



ff s (A)=2Tr 



P? s 



(Tr pf s /3f s j 



(27) 



Applying this to pi s as produced by the phase-flip chan- 
nel and using po given by Eq. (22) yields 



H s (\) 



4 (1 - rl) (r 2 - r 2 ) 



(1 - 2\y (1 - r 2 ) + 4A(1 - A) (1 - rf) 



(28) 



Thus the optimal single qubit, single channel use quan- 
tum Fisher information is attained when r, = and is 



H s opt (A) — 



4r 2 



1-(1-2A) r 2 



(29) 



Thus, for m independent channel uses the optimal quan- 
tum Fisher information is 



Hopt ind(A) 



4r 2 m 



1 - (1 - 2A) V 



(30) 



One physical implementation which yields a classi- 
cal Fisher information equal to that of -ff S opt(A) com- 
mences with po = (l + r&y^j /2 and follows a single 

channel invocation with a projective measurement in 
the basis {\+y),\-y)} where \±y) = (|0) ± i|l)) /V2. A 
straightforward derivation using the probability distribu- 
tion for the measurement outcomes and Eq. (7) results 
in F(\) = H e opt (A). 

Eq. (30) gives the optimal quantum Fisher information 
in the independent channel use protocol when all qubits 
have initial polarization r. This provides the bench- 
mark against which correlated state protocol estimation 
schemes will be compared. 



IV. CORRELATED STATE PROTOCOL 
PARAMETER ESTIMATION 



The results of section II C show that for pure input 
states, any correlated state protocol procedure cannot 
improve on the independent channel use protocol quan- 
tum Fisher information of Eq. (24). However, for mixed 
input states, the result of Eq. (30) does not saturate the 
absolute bound of Eq. (20). Here it is conceivable that a 
correlated state protocol procedure could improve upon 
the quantum Fisher information of Eq. (30) given the 
same number of channel invocations. 

Henceforth suppose that r < 1 and consider the cor- 
related state protocol procedure as illustrated in Fig. 3. 
Specifically suppose that po = (I + r& y )/2 and U prcp 
is as illustrated in Fig. 4. Denote the density operator 
for the system after the preparatory gate by p pre p := 
f^preppWprep- Then (see appendix A) 



(N-l)/2 
Pprop = ^ ] P(x) 



(31) 
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FIG. 3. Single preparation unitary scheme for improved pa- 
rameter estimation. There are n single qubits. A joint 
preparatory unitary, t/p rep precedes a phase-flip operation ap- 
plied once to each of the uppermost m qubits. This is of the 
form of the general scheme of Fig. 1 where To = U prC p, T m = I 
and Ti, . . . T m -i being swap gates between pairs of qubits. 



H — 



— H — 



FIG. 4. The preparatory gate, U plep , consists of a collection of 
controlled- Z operations applied to each distinct pair of qubits 
followed by a Hadamard transformation applied to each qubit. 
The preparatory gate is symmetric under interchange of any 
qubits since a controlled- Z gate is symmetric between the pair 
of qubits on which it acts. 



exception, if x ^ x' then pr x \ and p/ x <) have orthog- 
onal supports. The only exceptions to this are that 
P(x) = P(N — x)i a f ac t which is straightforward to show. 
Mathematically, the preparatory gate has engineered a 
splitting of the Hilbert space for the entire system into 
mutually orthogonal two-dimensional subspaces. 

The preparation scheme is symmetric under inter- 
change of qubits and without loss of generality, consider 
the situation where the phase-flip operation is applied 
once to each of the qubits corresponding to the m least 
significant (or rightmost) digits in the computational ba- 
sis representation. Then 



(iV-l)/2 
x=0 



(34) 



where pf f x ) results from applying the phase-flip opera- 
tion to p^) . Here 



\x)(x\ h> \x)(x\ 
\N -x)(N-x\ i-> \N-x)(N -x\ 

\x) (N - x\ H> (1 - 2A) m \x)(N- x\ 
\N - x) (x\ i— > (1 — 2A) m \N-x) (x\ 



(35) 



and thus 



Pi (a 



fix) + f(N - x) 



[x){x\+\N - x){N -x[ 



f(x)-f(N-x) 



+ - 2X) r 



\x){N-x\-\N-x)(x\ 



(36) 



where N = 2 n — 1 and 
-f(x) + f(N-x) 



P(x) 



+ i 



f(x)-f(N-x) 



\x)(x\+\N -x)(N -x\ 

\x){N - x\-\N - x){x\ 
(32) 



In Eq. (32), \x) = \x n ...xi) where x n ...x\ is the bi- 
nary representation of x, with the rightmost bit the least 
significant digit. Note that \N — x) can be obtained by 
flipping each bit value in \x). Finally 



(1 + r) j (1 - r) n ~ j 



(33) 



where j is the number of zeroes in the bit string for 
x. The support of pi x \ is the two dimensional sub- 
space spanned by {\x), \N — x)} and, with one type of 



Each of pi m has a two dimensional support and these 
are clearly orthogonal for distinct values of x in the range 
0, . . . , (N - l)/2. Thus Proposition 2 gives that H(X) = 



(N-l)/2 



l^x=0 



#(x) where # (x) 



Tr( 



dpi 



'-L (x) ) and L {x) is 



the score operator for p f ^ 

Proposition 1 applies to pf r x *, since it has two dimen- 
sional support and Tr pf ^ = f(x) + f(N — x) ^ 0. Di- 
rect calculation shows that Tr (p { — (Tr p { r x \) 2 = 
if and only if [fix) - f{N - x)f (1 - 2A) 2m = 
[f(x) + f(N - x)f . The fact that f(x) > implies that 
[fix) + f (N — x)] 2 ^ [f(x) - f(N - x)] 2 . Furthermore 
(1 - 2A) 2m s$ 1. This implies Tr (p { {x) ) 2 - (TV p { {x) ) 2 = 
if and only if (1 - 2A) = ±1 (i.e. A = 0, 1) and 
[f{x) + fiN - x)] 2 = [fix) - fiN - x)} 2 (i.e. fix) = 
or fiN — x) = 0). The latter is only true when r = 1. By 
the assumption that r < 1, Tr (/}f ( x )) 2 — (Tr/3f (a:)) 2 7^ 
and Proposition 1 gives that Li x ) derives from Eq (13). 
This together with Eq. (11) yields 
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if (A) = 4to 2 (1 



(JV-l)/2 
2A) 2m-2 £ 



a;=0 



[/(x)^/(iV- a; )] 2 [/(x)+/(jV^x)] 
[/(a?) + /(AT - x)] 2 - (1 - 2A) 2 ™ [/(*) - /(TV - x)] 2 ' 



(37) 



r 



The term within the summation of Eq. (37) yields the 
same quantity when x is replaced by N — x. Thus the 
range of the sum can be extended to to TV provided that 
the result is halved. Also the term within the summation 
only depends on the number of zeroes in the bit string 
for x and different x's can give the same contribution. 
Let j denote a possible number of zeroes. For x ranging 
from to TV the number of terms whose bit strings have 
j zeroes is (™). For each of these terms f(x)±f(N — x) 

[ 1 + r ) j (l - r) n - j ± (1 + r)"" J '(l - r) j ]/2 n . Thus 



H(X) 



(1 - 2A) 



2m -2 



2" 

x - f n 



c 2 d 



4 



(1 - 2A) 



2m „2 1 



(38) 



where 



Cj = (l + r) j (l-r) n - j 

_( 1 + T .)n-i(! _ r )J 

dj = (l + r) J '(l -r) n ~ j 
+(l + r) n - j (l-r) j . 



(39a) 
(39b) 



Equations (38) and (39) form the central result from 
which the principal conclusions of this article follow. 

The estimation accuracies of this correlated state pro- 
tocol and the optimal independent channel use protocol 
can be compared using the quantum Fisher informations 
of Eqs. (38) and (30). The constraints are that the polar- 
izations, channel parameter values and number of chan- 
nel invocations are identical for each protocol respectively. 
Given this, the fractional gain provided by the correlated 
state protocol is 



G(A) 



H(X) 



H op t ind(A) 



(40) 



and when G(A) > 1, the correlated state protocol per- 
forms more accurately than the independent channel use 
protocol. Equations (30) and (38) give 



this term is irrelevant in the sum). Thus as x increases 
from to 1, the fraction (1/x - r 2 )/(l/x m - in- 
creases. It follows that for a fixed polarization the mini- 
mum gain is 



1 



6V, 



On+1^2 / j \ j 
j=0 ^ J 





n \ c 

-r- if ' 



(42) 



if m 7^ 1 



and in each case these are attained when A = 1/2. Simi- 
larly for a fixed polarization the maximum gain is 



G, 



2 n + l r 2 j^\j) d j 1-^' 



(43) 



and this is attained when A = or 1. 

The other extremes of the region of parameter space 
are r — > and r — > 1. In the first case, equivalent to 
r«l, Cj « 2r(2j — n) and dj rj 2. Eq. (41) then implies 
that 



G mn(l - 2A) 2m " 2 as r -> 0. 



(44) 



As r — > 1, Cj — > unless j ' = or n. Specifically Co = 
-Cn -2" and d = d n 2". Thus Eq. (41) implies 
that 



G^m(l-2A) 2 



t _ 2 1-(1-2A) 2 
1 - (1 - 2A) 2 ' 



as r -> 1. (45) 



Note that this excludes the cases where both A = or 1 
and r = 1 here as Proposition 1 will give an SLD and 
Fisher information different to that of Eq. (37). These 
cases are beyond consideration here. 

The preceding results indicate important differences 
between situations in which the phase-flip operation is 
invoked only once and those where it is invoked more 
than once. 



G(A) 



m(l - 2A) 2 ' 

2n+l r 2 



3=0 



1 - (1 - 2A) 2 7 



2„2 



jjdj 1-(1-2A) 



2m„2 



(41) 



For a fixed polarization ^ r < 1 the gain is a mono- 
tonically increasing function of (1 — 2A) 2 . This follows 
from the fact that the gain function is a sum of terms of 
the form (1/x - r 2 )/{l/x m - c 2 /d 2 ) with x = (1 - 2A) 2 . 
Then, as shown in appendix B, c 2 /d 2 ^ r 2 (the only ex- 

and 



ception occurs when j = n/2, in which case c 



V. ESTIMATION WITH A SINGLE CHANNEL 

USE 

If the phase-flip is invoked only once then the corre- 
lated state protocol yields the following gain 



G ( A )-^Ti72 Ew)f r=7i 



1 - (1 - 2A) 2 r 2 



-\jjdj l-(l-2A) 2 c 2 /d 2 ' 



(46) 



As shown in appendix B, with just one phase-flip oper- 
ation and for any n ^ 2, G(A) > 1 if < r < 1. This 
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shows that, for single phase-flip invocation, the corre- 
lated state protocol of section IV gives a more accurate 
estimate of the phase-flip parameter than the independent 
channel use protocol whenever the polarizations used by 
each are identical. Additionally Eqs. (45) and (44) im- 
ply that G — > 1 as r — > 1 and that G — > n as r — > 0. 
The striking observation is that for sufficiently small po- 
larizations and a single phase-flip invocation, the corre- 
lated state protocol increases the estimation accuracy by 
a factor of approximately n regardless of the value of the 
parameter. 

For polarization values that do not approach these lim- 
its, a relevant measure of success would be the minimum 
gain, which is, according to Eq. (42), 




and is larger than 1. 

Typical plots of the gain, which illustrate the features 
described above, are provided in Fig. 5. These indicate 



Gain for n = 10 




FIG. 5. Gain using one phase-flip invocation. Both plots are 
over the range 0.05 ^ A ^ 0.95 and ^ r 1. 

that for a fixed value of A, the gain does not necessarily 
increase or decrease with increasing value of r, although 
it appears that G m - m decreases as r increases. Also, de- 
pending on the polarization and parameter value, it is 
possible to attain gains in excess of n. For example, 
with n = 2 the extremes are G m - in = 2/(1 + r 2 ) and 
G max = 2(1 + r 2 )/(l — r 2 ). Clearly a gain of arbitrar- 
ily large magnitude can be attained provided that, for 
example, A and r arc sufficiently close to 1. 



VI. ESTIMATION WITH MULTIPLE CHANNEL 

USES 

In the more general case of the protocol of Figs. 3 
and 4, there are m phase-flip invocations on n qubits. 
Here Eq. (41) indicates that G(l/2) = and thus there 
will be ranges of phase-flip parameters for which the cor- 
related state protocol provides a gain of less than 1 and 
thus gives a lower accuracy than the independent channel 
use protocol. 

A relatively simple example has n — 2 and m = 2 and 
the gain is illustrated in Fig 6. This indicates regions 




FIG. 6. Gain using two phase-flip invocations on two qubits, 
plotted for 0.05 sC A < 0.95 and < r < 1. 

in the parameter space for which G(A) > 1 and other 
regions for which G(A) < 1. Comparing Figs. 5 and 6 for 
n = 2 indicate regions for values of A near to or 1 in 
which the gain using two phase-flip invocations is larger 
than that using one phase-flip invocation. 

Whenever m ^ 2, Eq. (42) precludes any advantages 
for the correlated state protocol over the entire range of 
A. Similarly Eq. (44) results in the same loss of advantage 
over the entire parameter range when r — > 0. However, 
Eq. (43) shows that as A — > or A — > 1 the gain ap- 
proaches m times that of the correlated state protocol 
which uses only one phase-flip invocation. Thus for such 
parameter values the correlated state protocol will be 
superior to the independent channel use protocol. These 
facts arc illustrated in Fig. 7. 

It is possible to attain a substantial gain. For example, 
if r <C 1, then Eq. (44) shows that G(A) ^ n provided 
that 



The right hand side of this decreases as m increases. Note 
that a symmetrical situation occurs by replacing A by 
1 - A. 

Of special interest are the cases where the phase-flip 
is applied to all qubits (i.e. m = n). Here, the resulting 
upper bounds on A such that G(A) n arc illustrated 
in Fig. 8. This situation appears in NMR implementa- 
tions of quantum information processing, where all of the 
qubits arc subjected to the same or similar phase-flip op- 
erations. In room-temperature solution-state NMR im- 
plementations [14], the polarization of a single qubit is 
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Gain for m = 5 




0.8 0.2 



Gain for rn = 3 




0.8 0.2 



Gain for m = 1 




FIG. 7. Gain using various numbers of phase-flip invocations 
on five qubits, plotted for 0.0005 sC A sC 0.9995 and < r < 1. 



typically r sa 10 , which satisfies the assumptions re- 
sulting in Fig. 8. Considering, A = 0.10, a typical value 
appearing in Fig. 8, Eq. (3) shows that A ^ 0.10 when 
t ^ 0.22T2. In typical NMR quantum information pro- 
cessing experiments [31], Ti w 1 s and thus by limiting 
the period of evolution of the system to less than about 
0.2 s, the gain attained in estimating A will exceed n. 
This is certainly feasible to attain. 



VII. ENTANGLEMENT AND DISCORD 
CONSIDERATIONS 



It is customary to inquire about specifically which 
quantum resources could be responsible for gains of the 
type provided by the correlated state protocol. In other 
instances where quantum systems offer gains for pa- 
rameter estimation, attention has focused on the pres- 



A 

0.16 
0.12 
0.08 
0.04 








IB 



FIG. 8. Maximum value of A vs n such that G(A) ^ n for 
cases where r < 1. 



ence of entanglement or other types of quantum correla- 
tions [8, 26]. In the context of mixed state phase esti- 
mation it was shown that entanglement was not neces- 
sary for improved accuracy but it appears that quantum 
discord, a more general form of quantum correlation, is 
present whenever there is a gain in accuracy [26]. 

The preparation scheme for Fig. 4 is capable of pro- 
ducing states for which the measurement outcomes are 
correlated and it is of interest to determine the extent 
to which various quantum correlations are present in our 
scheme. The simplest context in which to address this for 
the correlated state protocol is that involving two qubits. 
For n = 2 the gain is 



G(A) = 2m(l - 2A) 2m " 2 (1 + r 2 ) 
1 - (1 - 2A) 2 r 2 



(l+r 2 ) 2 -4r 2 (l-2A) 5 



(49) 



We aim to determine whether the presence of entangle- 
ment is necessary for G(A) > 1 and whether other types 
of quantum correlation may be responsible for the gain 
attained by the correlated state protocol. We offer a 
preliminary investigation of these issues, focusing mostly 
on the presence of these correlations. A more complete 
treatment, which may consider issues such as the depen- 
dence of discord on the purity of the state and the extent 
of classical correlations in these protocols, is beyond the 
scope of this article. 

The existence of entanglement or other forms of quan- 
tum correlation is inferred from the density operator for 
the system. After m channel invocations, Eq. (36) yields 
that, in the basis {|00), |01), |10), |11)}, 



/I 



Pi = 



2 




1 — ir/i 





1 - r 2 
2 









1 - r 2 
2 









1 



(50) 
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where p, := (1 — 2 A)™. Note that the system density op- 
erator prior to channel invocation is the same as that for 
which A = 0, or, cquivalcntly, fi = 1. 
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A. Entanglement and separability 

A necessary and sufficient condition for a two qubit 
state to be separable (non-entangled) is that the partial 
transpose of the density operator for the state is a posi- 
tive operator [32, 33]. The partial transpose of pf is 



Pf 



\ 







2 

-ir\x 






irp 
1 -r 2 
2 








l + r 2 



\ 



(51) 
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This will be a positive operator if and only if its eigenval- 
ues are all positive. The block diagonal form of this ma- 
trix renders the computation of the eigenvalues straight- 
forward and it can be shown that these are positive if 
and only if < r < ^(1 - 2A) 2m + 1-|(1 - 2A) m | . Thus 
there will be polarizations such that for all A and to, pt 
is separable. Yet the results of Sections V and VI indi- 
cate the correlated state protocol yields greater accuracy 
than the independent channel use protocol for arbitrarily 
small polarizations and when m = 1. Thus the gains of 
the correlated state protocol cannot be ascribed to entan- 
glement. 



B. Discord 

Quantum discord is an alternative and more general 
characterization of the correlations between two quantum 
systems [34, 35]. Discord requires the density operator 
for the bipartite system, p and the reduced density opera- 
tors for the two the subsystems p A and /5 B , each of which 
is attained by taking the partial trace over the comple- 
mentary system. The quantum mutual information is 
l(p) := S(p A ) + S(p B ) - S(p) where the von Neumann 
entropy is S(p) := - Tr(/31og 2 p). 

A possible alternative approach to defining the quan- 
tum mutual information supposes that a measurement, 
described by POVM operators {Ek} where fc indexes the 
possible measurement outcomes, has been performed on 
system B and then inquires about the remaining informa- 
tion in subsystem A. Outcome k occurs with probability 
Pk := Tr (Ekp) and the state of system A following this 
outcome is p^ : = Tr b ( pEk )/pk- If the systems were in a 
product state prior to this measurement, then p A = pA\k 
for all k. Thus a measure of the classical correlation be- 
tween the subsystems, conditional on this particular mea- 
surement, can be quantified by S(p A ) — J2k PkS(pA\k), 
where the sum is an average over all possible measure- 
ment outcomes. Let S(p\{Ek}) := J2kPkS(f>A\k) ■ The 
overall classical correlation between the subsystems is 
them an maximum over all possible measurements and 
is C(p) := S(p A ) — min S(p\{Ek}) where the minimiza- 
tion is done over all POVMs on the subsystem B. The 



quantum discord is defined as the difference between the 
quantum mutual information and the classical correla- 
tion, Q := X(p) — C(p). In general, Q(p) ^ and it 
is invariant under local unitary transformations on the 
subsystems [34, 35]. 

The requirement for minimization over all possible 
measurements means that computing the quantum dis- 
cord is generally very difficult. However, analytical ex- 
pressions exist for certain classes of two-qubit states [36, 
37]. In particular an expression exists for for "X-states," 
which have the form, 



(52) 




where Cj are all real and a\ = a x , . . . . Let 



A 
Ai 
A 2 
A 3 



= 1 - Ci - c 2 - c 3 
= 1 - Cl + c 2 + c 3 
= 1 + ci — c 2 + c 3 
= 1 + ci + c 2 - c 3 



and 



Then 



Q{p) 



max{|ci| , |c 2 | , |c 3 |} 



1 3 1-c 

- 22 Xj log 2 Xj — log 2 (1 



(53a) 
(53b) 
(53c) 
(53d) 

(54) 



l + c 



log 2 (l + c). 



(55) 



The final density operator of Eq. (50) can be brought 
into the form (as it is the relevant values of Cj would not 
all be real) of Eq. (52) by the single qubit unitary 



U = 











(56) 



applied to each qubit. This will not alter the quantum 
discord. After this, direct calculation gives 

(57) 
(58) 
(59) 
(60) 

and 



(61) 



The discord for p pTep can be determined by setting A = 
in these expressions. Then c = r and 



A 


:= 1 - 


r 2 


Ai 


:= 14 


2r(l - 2A) m + r 2 


A 2 


:= 1 - 


2r(l - 2A) m + r 2 


A 3 


:= 1 - 


r 2 


c := 


= max 


{r 2 ,r(l-2A) m }. 



203) = 



1 



log 2 (1 



1 - r 



log 2 (l-r). (62) 



If r > 0, then ^2 > o since Q — when r = 0, 
this implies that the system state prior to the channel 
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invocation has non-zero quantum discord. In this sense, 
there are always quantum correlations present prior to 
channel invocation. 

After channel invocation, the discord depends on 
whether r > \i := (1 — 2A) m or not. In all cases, it is 
clear that replacing \i with —\i does not alter Q and thus 
it suffices to consider fi > (negative values are possible 
if m = 1). In appendix C it is shown that for fj, > 



— — > and — > 0. 

Ofi or 

With regard to A, one extreme is A = 
Here c = r and 

Q(p) = o. 



(63) 

1/2 (i.e. fj, = 0). 



(64) 



Thus regardless of the number of channel invocations 
there is no quantum discord when A = 1/2 (the final 
state of the system is clearly separable). However, with a 
single channel invocation, the gain in quantum Fisher in- 
formation is always at least 1 and approaches 2 as r —> 0. 
Thus the gains provided by the correlated state protocol 
cannot always be attributed to non-zero quantum discord. 

Beyond this extreme, i.e. fi > 0, the discord is non- 
zero. For a fixed value of A ^ 0, the discord is a monoton- 
ically increasing function of r. However, Eq. (49) can be 
used to show that the gain is not monotonically increas- 
ing. This requires finding r such that ^ = 0. Explicit 
calculation reveals that this occurs when r = or 



1 + v ± 2 v /(l + z/)V i 



4v r ' 



(65) 



where v = ( 1 
such that It = 

or 



- 2A) 2 . Table I lists examples where r is 
0. 



m 


A 


r 


1 


0.95 


0.66 


1 


0.99 


0.83 


2 


0.95 


0.48 


2 


0.99 


0.76 



TABLE I. Values of r such that |2 

or 



as calculated using 



Eq. (65). 



This indicates that the gain is not a monotonically 
increasing function of r. Thus there will be regions where 
the discord of the pre-measurement state increases while 
the gain attained decreases and an increase in the discord 
corresponds to a decrease in the gain. This is also evident 
from the lower graph of Fig. 5. 



first protocol, the independent state protocol, proceeds 
in such a way that the states of the system qubits are al- 
ways non-entangled. This is akin to the classical strategy 
of repeating the estimation procedure independently mul- 
tiple times and offers the same quantitative advantages. 
The second protocol, uses states which are correlated via 
an entangling unitary prior to channel invocation. If the 
number of channel invocations is identical and the polar- 
izations are identical, we have shown that the correlated 
state protocol can improve on the accuracy of the param- 
eter estimation when compared to the independent state 
protocol. The extent of the improvement depends on the 
actual parameter value but it can be substantially larger 
than the typical gain attained in the classical strategy. 

Our analysis demonstrates that the correlated state 
protocol attains advantages when the state of the system 
prior to measurement is separable and, in some cases, 
where there is zero quantum discord. Also, we have 
shown that, in some situations, an increase in quantum 
discord is associated with a decrease in the gain in es- 
timation accuracy. The gains afforded by the correlated 
state protocol cannot be attributed entirely to entangle- 
ment or non-zero quantum discord and it is unclear which 
quantum resources are possibly responsible for the gains. 

The correlated state protocol may be of use in esti- 
mating the dephasing time, Ti of a spin in NMR as 
this is closely related to the phase-flip parameter, A. 
However, we have assumed that the phase-flip param- 
eter is identical for all qubits; this is not necessarily true 
in implementations. For example, NMR systems used 
in quantum information processing provide different de- 
phasing times amongst the system qubits (for examples 
from NMR quantum information processing see [38-42]). 
It remains to be seen what gains methods like those of 
this article will yield in such cases. 

We stress that the correlated state protocol that we 
analyzed uses one particular preparatory unitary, C/ pre p, 
and we have not determined whether this is optimal. 

Finally, the methods used here are applicable to other 
Pauli-type channels, i.e. where a z is replaced by o~ n with 
n corresponding to any direction, in Eq. (1). These chan- 
nels are all related by single qubit rotations and the latter 
can be absorbed into the various parameter independent 
unitaries and measurements used in the estimation pro- 
tocol. 



VIII. CONCLUSION 
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Appendix A: State after preparatory unitary 



The Hadamard gate of Fig. 4 acts on a single qubit as 



7i £ o C-D™ l«) 

Vi=0 



(Al) 



and the controlled- Z between bits j and k via 

\z n ...zx) -> (-l) z ^ |z„...z!). (A2) 

Thus 

1 (A3) 



L/prep - 2 „/ 



where the binary representation of z is z n ■ • • z\ and the 
sum is over all possible values of these binary repre- 
sentations of z and y. Also z ■ y := J~]j—-t zjyj and 
s := z\Zi + ziz 3 + . . . + z n _iz n , in which each distinct 
pair of non-equal subscripts appears exactly once. 

In order to determine /5 prep := U pre pPiU^ Icpl it will be 
convenient to rc-cxprcss p\ ~ po ® po ® ■ ■ ■ ® po hi terms of 
the single qubit basis states, |±) = (|0) ± \l))/y/2. Here 



po = ?I+X+I+(i-?)I-X-I 



(A4) 



where q = (1 + r)/2. A convenient representation of pi is 
in terms of a sum of mutually orthogonal operators, each 
of which has a two dimensional support. Thus 



(JV-l)/2 

Pi= ^ (*) 

z=0 



(A5) 



where straightforward algebra using Eq (A4) gives 

A ( X ) = f(x)\x)(x\+f(N - x)\W^)(N~^\ (A6) 

with / as defined in Eq. (33). Here \x) is such that if 
the binary representation of x, then \x) con- 
tains the binary string for x with "0" replaced by "+" 
and "1" replaced by "— ". For example with n = 3, 

10} = !+ + +> = |+)|+>|+>, IT) = |+ + ->, |2) = |+ - +>, 

etc,. ... In this example p\ ( ) has support spanned 

by {|+ + +), I )}, pi (!) has support spanned by 

{|+ + — ), I h)} and so on. 

The preparatory unitary results in 



(JV-l)/2 

Pprcp = ^ ] P(x) 
k=0 



(A7) 



where 



P(x) '■— t^prcpPi (a')^prep 

= /(a;)C/p re p|^)(x|^ rep 

+f(N ^ x)U prcp \W~^)(N~^\Ul rcp . (A8) 

Each term in the decomposition of Eq. (A7) is an opera- 
tor with two dimensional support and these supports are 



all mutually orthogonal. Specifically the support of p( x ) 
is {[/p rep |3;), [/p r ep|iV — x)} and the unitary nature of the 
preparatory gate results in the orthogonality of distinct 
supports. 

It remains to determine an expression for U pTep \x) in 
terms of the computational basis (i.e. with bits strings 
consisting of "0" and "1" terms rather than "+" or "— " 
terms). For a single qubit 



') = 7f E (-1)*'" \*j) (A9) 



and this gives 



N 



z i= 



2™/2 



Ew zi+ - +Zn (-ir*i«>- 



(A10) 



Combining Eqs. (A3) and (A10) gives 



1 



N 



U P re P \x) = — J2 (i)^"^ (-1) 
y,z=0 



y-z-\-x-z-\-s 



\y) (ah) 



where s := z\Z2 + ^123 + . . . + z„_iz„, in which each dis- 
tinct pair of non-equal subscripts appears exactly once. 
We shall show that 



^prepl-^) 



1 + i 



1 - i 



\N 



(A12) 



(recall that \x) contains the bit string for x in terms of 
"0" and "1" while \x) contains the bit string for x in 
terms of "+" and "-"). Substitution into Eq. (A8) then 
yields the result of Eq. (32). 

In order to show Eq (A12), consider (x\U prep \x) . Then 
Eq. (All) gives 



\x-z-\-x-z-\-s 



1 N 

= ^Y,w i+ - +zn (- i y ( ai3 ) 



since (— l) x z+x z = 1. The term in the sum, 
(i) 21 "' z " (— l) s depends only on the number, m, of 
z n which are equal to 1. Clearly 7 := (i) zi+ "' +Zn = 



"-1 



The only bits which contribute to s are those for 



which both 



= 1 and Zk = 1. Gi 



that 



of 



Z\, . . . z n are equal to 1, the number of such such pairs is 
(™) = m(m — l)/2. Depending on whether this is even 
or odd, the term j3 := ( — l) s gives ±1. The various possi- 
bilities for to and their contributions arc summarized in 
Table II. 

This shows that whenever to is even the contribution 
to the sum is /3y = 1 and whenever it is odd, /?7 = i. As 
z ranges through all possible values, the number of times 
that to is even is 2 n /2 (exactly half of all bit strings have 



14 



m 


P 


7 


Pi 


4k 


1 


1 


1 


4k + 1 


1 


i 


i 


4k + 2 


-1 


-1 


1 


4fc + 3 


-1 


—i 


i 



TABLE II. Contributions to the sum of Eq. (A13). Here 



even numbers of bits equal to 1). Thus there are 2 n /2 
contributions of 1 and 2 ra /2 contributions of i. Thus 



(x\U p 



1 



(A14) 



Now consider (N — x\U plcp \x). Note that each bit in 
(N — x | is the bit flip of the corresponding bit in (x\. 
Then in Eq. (All) 



(N-x\(-ir\y) = (-l) 

Thus 

1 N 

(N - x\U wep \x) = - £(- 



(A15) 



>Zl+...+2„ 



(-l) s (A16) 



(x \ U p rC P \ % , 



which shows that (N — x\U plcp \x) = ^(a 
(1 — i)/2. This determines the two coefficients that ap- 
pear in Eq. (A12). Since U prep \x) has unit modulus, 
there cannot be any other components. This demon- 
strates Eq. (A12) and proves Eq. (32). 



Appendix B: Bounds on G(A) 

The general form of the gain is given by Eq. (41), which 
contains the term c^/d 2 . This enters into various bounds 
on the gain and following is a key result. 

Proposition 3: Suppose that r < 1. If j ^ n/2 then 



-L > r 2 
d 2 

with equality if and only if r — 0. If j = n/2 then 

c 2 

— = 

d 2 

Proof: Clearly if r = 0, then Cj = and thus Cj/dj 
Thus assume that r > 0. Let 



(Bl) 



(B2) 



f - ml 



(B3) 



By the definition of cj and dj, it follows that f n -j = ~fj- 
Thus it suffices to prove the result for j ^ n/2. 

Two special cases present themselves. First, if n is even 
and j = n/2, Eq (39) gives that cj — 0. Second, if n is 
odd and j = (n+l)/2, then let n = 21 + 1. Thus j = 1 + 1 
and direct substitution into Eq (39) gives Cj/dj = r. 

For all other cases j > (n + l)/2. Here 

dj(l - r) - 2(1 + r) n - j (l - r) j 



U 1 = 



d jr 



(B4) 



The numerator of this expression gives 



dj(l - r) - 2(1 + r)"- J (l - r) J = (1 + rY(l - r ) n - j+1 + (1 + r) n ~ j (1 - r) j+1 

-2(l + r) n - j (l-r) j 



(i + r y(i - r ) n - 3+1 - (i + r) n - j+1 (i - r y +1 . 



(B5) 



Now if < r < 1, then [(1 + r)/(l - r)] 2] 1 " > 1 pro- 
vided that j > (n+ 1)/2 (note that equality occurs if and 
only if r = 0.) Thus (1 + r) 2 ^ 1 "™ > (1 - r) 2 ^ 1 "". This 
implies that (l+rp(l-r) n ^+ 1 -(l+r)"^+ 1 (l-r)J+ 1 ^ 
0, with equality if and only if r — 0.. Furthermore, dj > 
and for r > 0, Eq. (B4) gives that fj — 1 > 0. This proves 
the proposition. □ 

It follows from proposition 3 that, if r < 1 then the 
summand in Eq. (41) satisfies 



c[ 1 - (1 - 2A) 2 r 2 
d~ 1 - (1 - 2A) 2m c 2 /d 2 



dj 1 



1- (l-2A) 2 r 2 
(1 - 2A) 2 "V 2 



with equality only possible if r = (note that for j - 
the result of proposition 3 is irrelevant since Cj 



(B6) 



n/2, 
= 0). 



Thus 



G(A) > 



m(l - 2A) 2m " 2 
2«+i r 2 

1 - (1 - 2A) 2 r 2 



1 - (1 - 2A) 



2m r 2 / ^ 
J=0 



In this context, the following is useful. 
Proposition 4-' If n ^ 2 then 



E 

J=0 



> 2 n+1 r 2 



-f. (B7) 



(B8) 



with equality if r = 0. 

Proof: If r = then Cj = and the equality holds. In 
general, straightforward algebra yields 

c 2 = d 2 ~ 4(1 -r 2 )™. (B9) 
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Thus 



n / x 2 " 



2™ +1 -4(l-r 2 )" V ( )t- (BIO) 



J=0 



n\ 1 



It will be shown that 



dj ^ 2(1 



„2\n-l 



(Bll) 



and thus 



4(1 -r 2 )" ]T 



n\ 1 



= 2 n+1 (l-r 2 ). (B12) 



This will establishes the proposition. 

It remains to prove Eq. (Bll). It is only necessary to 
show this for j ^ n/2 since d n -j = dj. First, for n = 2, 
there are two cases. With j = 1, d\ = 2(1 — r 2 ) and with 
j = 2, d 2 = 2(1 + r 2 ). Eq. (Bll) is clearly satisfied for 
both cases. Suppose that n > 2 and let 



(r) := 2 



(1-r 2 )' 



(B13) 

(r) = (1 - 
0. Suppose 

that j > (n + l)/2. Then we shall show that </>(r) is a 
monotonically decreasing function of r for < r < 1. 
First </>(0) = 1 and lim r ^i 0(r) = 0. Now consider 



Then if j = n/2, dj = 2(1 - r 2 )"/ 2 and 
r 2^n/2-i ^ ^ with equality if and only if r 



(1 -r 2 



4 



[r(2 - n)dj + (n - 2j) Cj ] . (B14) 



This is a continuous function of r and if <j)(r) is not mono- 
tonically decreasing or monotonically increasing, then ^ 
must equal zero at some point in the range < r < 1. 
Then = if and only if r{2 - n)d,- + (n- 2j)c,- . This is 
equivalent to r(2-n) (l + u 2 i-' n ) = (2j-n) (l - a 2 -?"™) 
where a = (1 — r)/(l + r). Clearly a < 1. But if n > 2 
and j > n/2 then the left hand side is strictly less than 
while the right is positive or 0. This is impossible and 
+ 1 0. Thus 7^ and has the same sign throughout 
the range < r < 1. So <f>(r) is either monotonically de- 
creasing or monotonically increasing. The latter is true 
since 0(1) < </>(0). This shows that <j>(r) < (f)(0) = 1 and 
proves Eq. (Bll). □ 



The results of proposition 4 together with Eq. (B7) 
imply that for < r < 1, 



G(A) > m(l - 2A) 



!m _ 2 1 - (1 - 2A) 2 r 2 
1 - (1 - 2A) 2m r 2 



(B15) 



It follows that with a single invocation of the phase- flip 
operation (i.e. m = 1), G(A) > 1 whenever < r < 1. 
Appendix C: Properties of Q 

Straightforward differentiation and algebra shows that 
9Q 



> 



(CI) 



whenever fi > 0. 
If 7' /j > then 

dQ r 

= - log 



dr 



(1 + r 2 + 2rfi)(l + r 2 - 2rfi) 



, V 1 
+ 2 !og 2 

+| log 2 
= 2 bg2 



(1 



(l + r 2)S 

+ 2r/j) 



(1 



2r/i) 



(1 + r 2 + 2r/x)(l + r 2 - 2r/z) 



(1 



(l + r 2)S 

+ 2r/j) 



(1 



(1 



-2r/i) 
2r/j) 2 



(1 + r 2 ) 2 



(C2) 



The argument of this logarithm is larger than or equal 
to 1 (only possible if r = or fi = 0). This shows that if 
r > fi then ^£ > 0. 

On the other hand if fi > and r ^ ft then 



as 

Sr 



7: lop 
2 e 



(1 + r 2 + 2r/i)(l + r 2 - 2r/i) 



^ 1 

"2 bg2 



(1 + r 2 ) 2 
(l + r 2 + 2r/i)(l-r M ) 



> 2 lo §2 



(1 + r 2 - 2r/j)(l + r/x) 
(l + r 2 + 2r/u) 2 (l-r/i) 



(C3) 



Let 



Then 



(l-r 2 ) 2 (l + r A t) 

(1 + r 2 + 2r/j) 2 (l - r/j) - (1 - r 2 ) 2 (l + rfi). (C4) 

u = 2r [/i(l - r 2 )(l + r 2 ) + 2r(l - r 2 /i 2 ) 

+2r 2 /j(l - mu 2 )] (C5) 

and for every ^ r, mu ^ 1 every term here is positive. 
Thus to > and the argument of the last logarithm is 
larger than 1. Thus if fi > and r ^ fx then ^ > 0. 
This shows that in all cases ^ > 0. 
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